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We report experimental observations that a bubble rising in water in the presence of a sound field is signifi- 
cantly slowed down, even at moderate acoustic pressures. 

We measure the mean rise velocity of bubbles under various acoustic forcings and show this effect occurs 
if the noise spectrum matches or overlaps bubble resonance. We render surface oscillations and translational 
movements of bubbles using high speed video imaging and thereby identify Faraday waves on the bubble wall 
as the cause for the velocity reduction. The associated mechanisms are discussed in terms of induced forces. 
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The velocity at which buoyant gas bubbles rise to the 
surface of a liquid is key to geophysical and industrial 
processes^. Hence, considerable effort has been devoted 
to understanding the hydrodynamics of a single bubble^, 
underlying the major role played by the shape and the 
wake on its buoyant rise. 

Bubbles are also known to be highly sound-sensitive 
objects^; a property which also has implications in vari- 
ous fields: from sound propagation in oceans^, to sono- 
luminescence and sonochemistry^, or even for biomedical 
applications^. Acoustically, a bubble responds to acous- 
tic pressure fluctuations as a peculiar resonator since its 
resonant Minnaert frequency /m (Hz) obeys 

fM w 3.26/i? , valid for bubble radius R > IQ-^ m, 

(1) 

and corresponds to a forcing wavelength in water much 
greater than the bubble size, such that the bubble un- 
dergoes radial oscillations. Yet, beyond a critical accel- 
eration of the bubble wall, the parametric Faraday insta- 
bility is triggered and shape oscillations (surface waves) 
developii^. Furthermore, because a bubble is a sharp os- 
cillator, this threshold is lowest in terms of forcing pres- 
sure if the bubble is excited near resonance. Acoustic 
forcing of these waves has interesting applications such 
as bubble sizing^, electrodeposition patterningiS or mass 
transfer enhancementsii. 

Although most of the above cited applications involve 
buoyant bubbles rising in noisy environments, fewer stud- 
iesi^"— concern these coupled aspects. The shape of a 
bubble being key to its hydrodynamics, surface waves are 
expected to significantly modify the bubble rise. How- 
ever, existing results focus on the influence of high inten- 
sity sound fields on bubble motion, in which case effects 
of the Faraday instability, if triggered, overlap with other 
mechanisms such as the acoustic radiation force. 

Hence, it is logical to address the following question 
: "what is the effect of a low intensity, permanent and 
resonant sound field on the rise of a bubble?" 

In this letter, we report experiments studying the in- 
teraction of a single rising bubble with a resonant sound 
field. Rise velocities with and without sound are com- 
pared by time of flight measurements to capture effects 



of the acoustics. Our results are discussed to identify the 
physics at play but the complete theoretical modeling is 
beyond the scope of this letter. 
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FIG. 1. Right : Sketch of the experimental set-up featuring 
the bubble-on-demand system, the 4 lasers that detect bubble 
crossings to measure rise velocities and to trigger events, the 
sound emission chain and the digital cameras. Left : the time 
sequence of a bubble rise is decomposed in key steps including 
injection, sizing, velocity measurements and sound emission 
triggering. 

The experimental set-up is sketched on Fig. [TJ Air 
bubbles are produced at the bottom of a tall Plexiglas 
water tank: a cubic visualization cell (9 cm wide) linked 
to a round tube (6 cm i.d, 40 cm height) leads to a 
squared vessel (30 x 30 x 90 cm). Experiments are carried 
out using purified deionized water at room temperature 
(T « 25" C). The vessel is sealed to limit contamination. 
Single bubbles are released from a capillary silica tube 
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(20 i.d., 90 fj,m. o.d.) connected to a fixed pressure 
source. A home made system (by-pass and electrovalve) 
enables on-demand bubble production every 90 s to en- 
sure that the previous bubble and its wake have van- 
ished. Upon release, the bubble is detected and sized 
at the bottom of the tank by a digital camera synchro- 
nized to a stroboscopic flash. The bottom bubble radius 
Rb is extracted by live image processing using a Mat- 
Lab routine (1% accuracy). Very good reproducibility 
is obtained since the diameters of successive bubbles are 
identical within the measurements error margin. 

During rise, 4 lasers (3 lines and 1 sheet) placed at 
heights 5 cm, 53 cm, 93 cm and 129 cm above the injec- 
tion nozzle detect bubble crossings, giving 4 trigger times 
and thus 3 rise velocity measurements by time of flight 
(Vf,, Vm and Vt for bottom, middle and top respectively - 
see Fig. [1] for details). Crossing detections are also used 
as triggers : for sound emission and high speed camera 
acquisition by the laser 3 and for sound deactivation by 
the top laser sheet. 

Sound is produced by an ITC 1001 spherical transducer 
(an omnidirectional projector) powered by an audio am- 
plifier (HPA A2400) followed by a transformer (ratio 4.5). 
Monochromatic or band filtered white noise signals can 
be generated using a B&K 3310 Pulse system. The pro- 
jector is placed at the bottom of the top vessel in most 
experiments. Sound pressure signals are recorded by a 
hydrophone (B&K type 8103) connected to a charge am- 
plifier ( B&K Nexus). The tank being a resonant cav- 
ity, the sound field along the vertical centerline of the 
tank is mapped prior to each series of experiments in- 
volving monochromatic sound. The hydrophone is then 
placed near a wave crest 2 cm away from the center- 
line to give a measure of the sound pressure in the res- 
onant tank throughout the experiment^. A high speed 
video camera (Photron SAl) captures the rise and os- 
cillations of the bubble in synchronization with the hy- 
drophone sound pressure signal using a dedicated syn- 
chronized fast- acquisition system (Photron MCDL). 

The experimental protocol is the following: First, bub- 
bles are produced until stability of the injection nozzle is 
reached, setting a reference bubble radius for the rest 
of the series. Assuming an isothermal hydrostatic expan- 
sion of the bubbles, the radius Rt at the height of laser 3 
is predicted and injected into a refined version of ([T]) {ci^ 
p. 183.) to calculate the resonant frequency /m, taking 
dissipative effects into account. Reverberation effects can 
be shown to be negligible in our configuration^^. After- 
wards, a series of bubbles are released without acoustic 
forcing to determine the three reference quiescent rise ve- 
locities {V^, and V^). Only then is sound emission 
activated and the effect of acoustics are probed. The 
bottom Vb and middle Vm velocities are still measured in 
a quiescent environment. Along with the bubble radius 
Rb, they are indicators of reproducibility as they should 
remain equal to the reference values i?" , and V^^ , and 
are used to rule out random events such as changes in 
bubbling regime. Finally, the effect of acoustics on the 



bubble ascent is captured by the top rise velocity Vj, 
measured above the laser 3 that triggers sound emission. 

We run experiments in which bubbles are forced at 
a fixed frequency / near resonance /a/ and at varying 
sound levels. The top rise velocity Vt normalized by the 
reference value is plotted in Fig. [2] as a function of the 
acoustic pressure P for both monochromatic sound and 
band filtered white noise of fixed frequency span. The 
curves featured on Fig. [2] are typical of all conducted ex- 
periments, in which bubble radii ranged from 300 fj.m to 
500 fim. Firstly, we observe that the acoustic forcing at 
resonance reduces rise velocity by up to 30% for moder- 
ate noise levels of order P sa 100 Pa, and by more than 
50% for P > 300 Pa. Secondly, there is always a clear 
pressure threshold Pc upon which acoustic effects come 
into play, around P = 50 Pa for the plotted measure- 
ments. This threshold does not differ significantly from 
monochromatic to band filtered white-noise cases. Fur- 
thermore, we recover the shape of the curves whatever 
the projector's position, indicating that only the value 
for the threshold Pc varies with the acoustic forcing char- 
acteristics (projector position and spectrum). 




FIG. 2. Bubble velocity Vt normalized by the quiescent ref- 
erence velocity V^, as a function of the root mean square 
(RMS) sound pressure P for monochromatic (+) and band 
filtered white noise (o) emission. For the monochromatic case 
(+), R'i = 452 fim, Vt° = 19.2 cm-s-\ /a/ = 7.04 kHz and 
/ = 7.2 kHz. For the filtered white noise, = 446 fj,m, V^' = 
18.9 cm-s~^, /a/ = 7.13 kHz, central frequency / = 7.2 kHz 
with a span A/ = 1.6 kHz. Pressure bands relate the differ- 
ent surface oscillation regimes to the velocity decrease. 
Insert: Vt/Vt' as a function of the acoustic frequency / for 
a monochromatic forcing: = 447 fim, = 19.0 cm-s^^, 
/a/ = 7.12 kHz, P = 114 Pa. 

To test the frequency sensitivity of the phenomenon, 
we excite a series of bubbles by monochromatic forcings 
at a constant sound pressure P above the threshold Pc for 
velocity reduction and at various frequencies around reso- 
nance. Insert in Fig.[2]shows that bubbles arc only slowed 
down by acoustic forcings near resonance for a narrow 
bandwidth of 300 Hz approximately, which is of the order 
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of the bubble's resonant bandwidth A/o = 280 Ha^. We 
also verify that forcing at harmonics and subharmonics 
of the resonant frequency have no effect on the rise veloc- 
ity. We obtain similar results using band-filtered white 
noise forcings but the added parameter of frequency span 
leads to a more complex parametric analysis. 

Hence, the velocity reduction induced by the acous- 
tic forcing is shown to be a threshold effect restricted 
to a narrow frequency band around resonance. Other 
phenomenons involving acoustic forcing of bubbles share 
such onset properties, such as subharmonic emissions of 
resonant bubbles^ and mass transfer enhanccmentsii. In 
these problems, the origin for the threshold behavior is 
the parametric Faraday instability, which triggers sur- 
face waves on the bubble wall. Theoretical models and 
experimental measurements^iii predict the triggering of 
Faraday waves on bubbles of 500 fim radius for sound 
pressures of approximately P w 5 • 10~^/i? = 100 Pa, 
which is of the order of the threshold Pc we measure for 
rise velocity reduction. To substantiate the link, we vi- 
sualize in detail bubble behavior under acoustic forcing 
at resonance using high speed video imaging (62500 fps). 

In the case of monochromatic sound emission, the 
videos reveal a clear transition in bubble surface motion 
at the pressure threshold Pc- Below this onset value, bub- 
bles oscillate radially while above, they undergo shape 
oscillations. The bubble geometric distortions are first 
localized near the equator, as shown in figure [3l^a). Such 
equatorial modes have been observed for tethered bub- 
bles and modeled in the framework of Faraday waves as 
self-focusing ripplesi^. Here, this particular feature is 
visualized on rising bubbles. 
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FIG. 3. Frame sequence for rising bubbles in (a) equatorial 
mode and (b) mode 8 oscillations. Curves below represent the 
instantaneous acoustic pressure p{t) near the bubble ( - ) and 
the normalized projected area fluctuations S' of the bubble 
(o). Rb = 420 f = hi = 7.6 kHz, (a) P = = 15 Pa 
and (b) P = 31 Pa, Vt° = 17.8 cm-s"\ frame rate 62500 fps. 

As the acoustic pressure increases to a value of P w 
2 • Pc, the pattern formed by the Faraday waves changes 
from equatorial distortions to 8-mode oscillations, as 
shown on figure [3l[b). This observation is consistent with 
theoretical results^iiii, which predict the 8-mode as the 
most unstable parametric surface wave for a non moving 



oscillating bubble of size R = 420 ^im. 

Well above the pressure threshold, for 4 < P/Pc < 6 
approximately, multiple unstable modes interact and give 
rise to less regular geometric deformations and complex 
patterns. For even stronger forcing of typically P/Pc > 6, 
the bubble oscillations are highly distorted, no longer ax- 
isymctric nor harmonic. In this regime, the rise velocity 
Vt decreases to 50% of the reference value and the 
bubble "dances": it follows an erratic path along the 
vertical and fragmentation occurs, as often observed in 
acoustic trapping experiments^^. Each regime is illus- 
trated by a movie in the auxiliary materiali^. 

In the video sequences under white noise emissions 
band-filtered around resonance, the bubble exhibits non- 
stationary volume and shape oscillations, as both the 
local frequency and amplitude of the sound wave that 
impinges on the bubble are time dependent. Therefore, 
the onset threshold for Faraday waves is only met tran- 
siently and the instability never reaches a steady stateSS. 
Nonetheless, the surface oscillations of the bubbles can 
still be qualitatively described as sequences of equatorial 
ripples, n-mode oscillations and mode superpositions. 

Additional quantitative evidence is brought by the fre- 
quency of the surface oscillations which can be estimated 
from the video sequences and the time traces of the bub- 
ble's projected area and compared to the sound pressure 
signals, shown in Fig. [31 It is clearly subharmonic of the 
forcing frequency, a key feature of the Faraday instability. 

As a whole, these observations establish surface Fara- 
day waves as the cause for the rise velocity reduction 
but the underlying mechanism needs to be discussed. In 
our experiments, the bubble is sufficiently small that, in 
the absence of sound, it rises along a vertical lino^ at 
a velocity fixed by two competing forces: the buoyant 
driving force Fh = pVg in which p is the density of wa- 
ter, V the bubble volume and g gravity, and a drag force 
Fd = — l/2Crf/97rP^ \V\ V where Cd is a shape-dependent 
drag coefficient^! and V the bubble velocity. 

When sound is produced below the Faraday threshold 
Pc, the bubble undergoes radial volume oscillations (a few 
percent in amplitude: SR/R cx 10~^), leading the force 
balance into a non-steady state. Hence, the rise velocity 
also oscillates around a mean value, such that the un- 
steady added mass^^ and history^^ forces come into play. 
Yet, we observe that the mean rise velocity is not altered. 
Indeed, the contribution of these forces to the mean mo- 
mentum balance is at most 0{5R? / R^) when compared 
to the buoyancy force. Nevertheless, other effects have 
to be examined: First, velocity fluctuations also couple 
with radial oscillations, giving rise to a streaming effect. 
Under the operating configuration, the streaming flow 
is of order— 1 mm-s~^, which is negligible compared to 
the bubble rise velocity . Second, the sound field pro- 
duces a radiation pressure on the bubble^^. However, this 
force requires strong acoustics fields to be significant, as 
illustrated by the following scaling argument. Pressure 
fiuctuations of amplitude P' induce volume oscillations 
V cx R/{p£,P) ■ P', where ^ cx 10^^ is the damping co- 
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efficient of the bubble and / tfic forcing frequencj*^. Tlie 
pressure gradient scales as VP oc f /c ■ P' with c the 
sound velocity in water. At resonance, the forcing fre- 
quency / is given by ([1]), such that the radiation force 
relative scales as Fr/Fb oc< V • VP > /F^ oc IO'^PVP 
relatively to buoyancy. Hence, the radiation force is sig- 
nificant for acoustic pressures above Pr ~ IO^/a/P be- 
yond the threshold for Faraday waves Pc ~ 5 • 10~^/P, 
provided bubbles are at least micron sized. 

For sound pressures above Pc, the same forces com- 
pete but surface waves now superimpose onto the volume 
oscillations. The exact mechanism by which the Fara- 
day waves reduce the rise velocity is difficult to identify 
given the number of forces acting upon the bubble and 
the complexity of the shape oscillations. A likely candi- 
date is the increase in the mean hydrodynamic drag Fd 
through modification of the drag coefficient Cd- Such 
variations are be significant: the drag coefficient of a 
cube can be twice that of an equivalent sphere, depend- 
ing on its orientation. Here, given the observed angular 
shapes, Cd always increases when the bubble departs its 
spherical geometry, such that the mean drag coefficient 
also increases. Nevertheless, other mechanisms cannot be 
ruled out, since, for example, Faraday waves have been 
shown to greatly enhance acoustic streamings^. Hence, 
new models and numerical simulations are necessary to 
the full understanding of the observed rise velocity re- 
duction. 

To conclude, our observations clearly establish the 
strong impact of surface Faraday waves on a bubble's 
buoyant rise. Such waves are most easily triggered by 
resonant acoustic forcings such that moderate sound lev- 
els can induce important rise velocity modifications. This 
suggests to take into account the acoustic noise in multi- 
phase flow problematics: air entrainment in oceanic surf 
zones^l, charge losses in bubbly flows, or turbulence drag 
reduction by microbubble injection^^. From a more ap- 
plied point of view, our observations could initiate new 
strategies for industrial processes where the bubble rise 
rate is crucial, or for soft-matter sorting by acoustophore- 
sis. Indeed, as shown in a recent paper—, some cancer 
cells have a lower membrane elasticity and this contrast 
in stiffness should modify surface waves triggered on the 
cell wall, which in turn modifies the flow properties of 
the cell. It is the subject of ongoing investigations. 
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